The role of initial entanglement and nonGaussianity in the decoherence of 
photon-number entangled states evolving in a noisy channel 
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We address the degradation of continuous variable (CV) entanglement in a noisy channel focusing 
on the set of photon-number entangled states. We exploit several separability criteria and compare 
the resulting separation times with the value of non-Gaussianity at any time, thus showing that in 
the low-temperature regime: i) non-Gaussianity is a bound for the relative entropy of entanglement 
and ii) Simon' criterion provides a reliable estimate of the separation time also for nonGaussian 
states. We provide several evidences supporting the conjecture that Gaussian entanglement is the 
most robust against noise, i.e. it survives longer than nonGaussian one, and that this may be a 
general feature for CV systems in Markovian channels. 
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Continuous variable (CV) quantum information has 
been developed with Gaussian states and operations 0- 
3]. However, in the recent years also the non-Gaussian 
sector of the Hilbert space has been taken into consider- 
ation. This interest is due to the potential role of non- 
Gaussianity in enhancing long-distance quantum com- 
munication based on entanglement distillation 0, [B| and 
swapping, quantum memories Q , cloning Q and telepor- 
tation |8| . In turn, it has become of interest to analyze 
non-Gaussian states in realistic conditions [jij], where de- 
coherence due to dissipation and thermal noise unavoid- 
ably leads to degradation of entanglement. Our work 
is indeed motivated by the following general question: 
in case of transmission through a noisy environment is 
there any advantage in using non-Gaussian states? Do 
they lose entanglement in a longer time? We provide 
evidence for the answer to be negative, thus supporting 
the conjecture that Gaussian entanglement is extremal 
in terms of robustness against decoherence due to noise 
and dissipation. 

In order to address the above questions, in this Letter 
we consider a broad and meaningful class of CV bipartite 
states endowed with perfect correlations in the number of 
photons: photon-number entangled states (PNES). The 
latter have Schmidt decomposition in the Fock basis, i.e. 
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with real coefficients tp n £R,ip n > 0, Yl^Lo = The 
advantages of considering these states are twofold. They 
are sufficiently simple for analytical study, and at the 
same time meaningful since several experimental realiza- 
tions have been reported and quantum communica- 
tion schemes involving PNES have been proposed 11 1. 



Furthermore, the set of PNES contains mostly non- 
Gaussian states but includes (as a subclass) two-mode 
squeezed vacua, i.e. the basic Gaussian resource for CV 



quantum information, thus allowing for a direct compar- 
ison between Gaussian and non-Gaussian states. Finally, 
PNES are good candidates for long-distance quantum 
communication, because they have been already proved 
robust against some kind of noise, e.g. phase diffusion 
[l2| . We consider several special subclasses of PNES with 
specific parametric dependence, as well as randomly- 
generated [l3| (truncated) PNES, in order to draw some 
general conclusions about the typical behaviour of en- 
tanglement dynamics. In particular, we focus on ran- 
dom PNES with decreasing profile (i.e., tp n > ip n +i) 
and on the following parametric subclasses (we omit nor- 
malization): (i) the two- mode squeezed vacua or twin- 
beam states (TWB) ip n cx x n < x < 1 which are 
the sole Gaussian states within the PNES class and rep- 
resent the preferred (Gaussian) resources in protocols 
involvin g C V entanglement; (ii) the photon subtracted 
(PSSV) [HI ip n cx (n+l)x n+1 and the photon-added two- 
mode squeezed vacua (PASV) fHf ip n cx to" -1 , which 
are obtained from the TWB by the experimentally fea- 
sible operations of photon subtraction g — > a\a2Qa\a\ 
and addition g — > a\a\gaia,2 respectively [1(3]; (iii) the 
pair-coherent or two-mode coherently correlated states 
(TMC) with Poissonian profile t/>„ 

mean energy of PNES is (ii\a{ai 
N = X)^Lol' ( /'» 1 l 2n > whereas correlations between the 
modes can be quantified by C = Re^^°_ ip^ tp n +i (n+ 1) 
and entanglement is given by the Von-Neumann entropy 
of the partial traces eo = — J2 n 1°6 ■ ^ n turn, the 
covariance matrix (CM) of a PNES equals that of a sym- 
metric Gaussian state in standard form, with diagonal 
elements equal to N + \ and off-diagonal blocks given by 
C = diag(C,-C). 

The propagation in noisy channels can be modelled as 
the interaction of the two modes with two independent 
thermal baths of oscillators. The resulting dynamics is 
a Gaussian channel, governed by the two-mode Master 
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equation (ME) 

e = E ^ ^K 1 ^ + ?(^' + !) ( 2 ) 

describing losses and thermal hopping in presence of (lo- 
cal) non-classical fluctuations of the environment. Dot 
stands for time-derivative and the Lindblad superoper- 
ator is defined by L[0]g = 20gO^ - O^Og - gO^O. 
r is a loss coefficient and Nj are the mean photon- 
numbers in the stationary state, which is a thermal state. 
We consider baths at equal temperature iVi = JVo = 
Nt- The above ME admits the operator solution flij : 
g(t) = A t g{0) = Tr 34 [f/ t (f>(0) ® v 34 )U t ], where A t de- 
notes the evolution map corresponding to the noisy chan- 
nel; 3, 4 are two additional fictitious modes in a ther- 
mal state ^34 = v 3 ® v 4 \ U t = Ui 3 (( t ) ® ^24 (Ct) an d 
Uij((t) — Gxp(Qa\aj — Qa^ai) is the two mode mixing 
operator, with £ t = arctan(e rt — l) 1 / 2 . Using this so- 
lution, the evolved density matrix g t can be computed 
numerically from the initial state go upon truncating the 
Hilbert space dimension. In our study we consider states 
with total energy < 2N < 10 and dimension D = 20. 
In this range of energies, and for all subclasses of states, 
the adopted truncation results in a negligible error. We 
emphasize that the map A t , being the product of two lo- 
cal maps, can only disrupt quantum correlations: for any 
Nt ^ we have a complete loss of entanglement within 
a finite, state dependent, time t s = t s (g) which we refer 
to as the separation time. 

In order to estimate t s for non-Gaussian states sub- 
jected to the action of A t we make use of several en- 
tanglement criteria and this also enables a comparison 
of their performances in detecting entanglement. As it 
is well known, in the CV case a necessary-and-sufncient 
separability criterion exists only for Gaussian states [lj| : 
Simon's criterion (SI) for separability is equivalent to the 
positivity of the partial transpose density matrix and 
says that a Gaussian state is separable iff d_ < 1/2, 
where d_ is the least symplectic eigenvalue of the CM 
of the partial-transposed state. When dealing with CV 
non-Gaussian states, Simon's criterion (which is equiv- 
alent to the separability of a Gaussian state having the 
same CM as the given state) is only sufficient for en- 
tanglement. This actually holds for any available cri- 
terion: if the state is entangled, a given test may or 
not detect its entanglement; in turn, if no test detects 
entanglement, we can not conclude separability of the 
state. The Simon separation time t SI can be computed 
analytically. At the level of CM, the map A t induces 
the evolution cr t = aoe~ rt + ctoo(1 — e _rt ), where a oo = 
diag(A r + 1/2, . . . , N T + 1/2) is the asymptotic thermal 
state CM. The CM of the partial-transposed state is given 
by Acr t A, where A = diag(l, 1, 1, —1) [lit ], and we have 
d- = (N T + 1/2)(1 - e- rt ) + (N + l/2)e- rt - |C| e - rt . 
Therefore, for N T — PNES are entangled at any time, 



whereas for N T ^ Owe have a lower bound to separability 

Besides Simon's criterion, we will make use of three 
different criteria which provide independent separabil- 
ity conditions. The first is the extension of SI given 
by Shchukin and Vogel [20] (SH) based on the evalua- 
tion of a series of M x M matrices whose entries are 
moments up to a given order: non-positivc-dcfinitencss 
of any finite submatrix is a sufficient condition for en- 
tanglement. By considering the minor defined by the 
first and second-order moments only (M — 5) we ob- 
tain a condition which is equivalent to SI. If we con- 
sider larger minors, moments of higher order are involved 
and we get a stronger condition. Here we consider mo- 
ments up to order 8. The second criterion has been in- 
troduced by Sperling and Vogel (SP) [2l| and it is based 
on linear entanglement witnesses. A state g is entan- 
gled if (4>\g\4>) > max n {|m ra | 2 } where \<j>) is a pure entan- 
gled state with Schmidt coefficients {m„}. We test this 
condition by using 10 4 randomly generated witnesses of 
the form \<f>) = X^ =1 (f> n \n)\n) with D = 20, i.e. the 
witnesses are themselves truncated random PNES. This 
form is chosen since the bath does not create quantum 
correlations but only destroys those originally present. 
Finally, the realignment criterion [53] (RE) is based on 
positivity of a linear contraction map: a state g is entan- 
gled if ||g|| > 1 where \ \A\ \ denotes the trace norm of op- 
erator A and (m\(^\g\n)\v) = (m\(n\g\v)\^t) . Using these 
criteria, we obtain lower bounds on separation times. In- 
deed, for any given criterion K and state g, let us denote 
by t K (g) the maximum time for which K proves that g is 
entangled: clearly t K (g) is a lower bound for t s . Consid- 
ering the best bound we have t s {g) > t M = m&x K t K (g). 

The propagation in noisy channels, besides entangle- 
ment, also destroys the non-Gaussian character of the ini- 
tial state, which unavoidably evolves towards the asymp- 
totic, Gaussian thermal state. We shall take into ac- 
count both processes (separation and Gaussification) in 
parallel and explore the relations between them. The 
non-Gaussian character of a state g is measured by 
6(g) = S (r) — S (g) i.c, the relative entropy between g and 
the reference Gaussian state r having its same covariancc 
matrix [23[. In order to explore the effect of noise in a 
wide range of conditions and initial states we consider 
TWB, PSSV, PASV, TMC and random PNES of differ- 
ent energies and compute the evolved density matrix for 
< t < 15 in units of inverse loss 1/r. At any time t, 
entanglement is tested with all the above mentioned cri- 
teria and the value of the non-Gaussianity 6 is computed. 
From these data we evaluate t K , i.e. lower bounds to sep- 
aration times according to different criteria, and Gaussi- 
fication times t G , i.e. times for which non-Gaussianity 5 
falls below a fixed Gaussification threshold S G (we con- 
sider different thresholds S G = 0.1,0.01,0.001). The pro- 
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cedure is then repeated for different values of the temper- 
ature T corresponding to N T in the range [10~ 5 , 1CP 1 ]. 
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FIG. 1: (Color online) Separation and Gaussification times 
for TMC as a function of the mean energy for low (N T = 10~ 5 , 
left) and high (N T = 10 _1 , right) temperature. In both plots 
we report separation times according to different criteria: taB 
(green, circle), t s p (blue, square), t S n (red,triangle), t S i (pur- 
ple,rhombus), and Gaussification times for different thresh- 
olds: S G = 10 _1 (solid black, triangle), S G ~ 1CF 2 (dashed 
black,triangle), 5 G = 10~ 3 (dotted black,triangle) . 



We start describing the results of our analysis by fo- 
cusing on TMC. In Fig. Q] we report t K for TMC and 
different criteria as a function of N for the lowest (high- 
est) temperature considered N T = 10~ 5 (10^ 1 ). It turns 
out that at any temperature SI, SH and RE criteria yield 
similar curves whereas the SP criterion works only at low 
T. We point out two general features: (i) t K is a decreas- 
ing function of T, i.e. entanglement is strongly corrupted 
as the temperature increases; (ii) both at high and low 
T, t K rapidly increases to an asymptotic value t K which 
is reached at N ~ 1/2 and then remains almost constant. 
In Fig. [1] we also show Gaussification times of TMC as 
a function of energy. We see that the behaviour of non- 
Gaussianity is only weakly affected by the increase of T. 
Upon comparing separation and Gaussification times we 
notice that at low T states become nearly Gaussian well 
before they become separable: t G < t K < t s . At high T, 
on the contrary, Gaussification times are greater than our 
bounds on separation times: t G > t K . The analysis of all 
other PNES subclasses (TWB, PASV, PSSV and random 
PNES) reveals the same qualitative behaviour described 
for TMC for the evolution of both entanglement and non- 
Gaussianity (with the obvious exception of TWB, whose 
non-Gaussianity is always zero). In summary, we have 
numerically proved that for the whole class of states we 
have considered and at any temperature SI, SH and RE 
criteria yield qualitatively the same results. In addition, 
Simon's criterion, which offers analytical advantages, is 
the optimal one. Furthermore, the separation time de- 
creases with T, and the dependence on the energy can be 
appreciated only for small N, while they quickly reach 
their asymptotic values as N increases. As for the non- 
Gaussianity, we have t G < t K < t s at low T and t G > t K 



at high T. This deserves futher consideration. Indeed, 
the relation t G < t K < t s suggest that for low T the 
bounds provided by Simon's criterion properly estimate 
the actual PNES separation times i.e., t SI ~ t s . This 
can be understood by first noticing that when t > t G 
the states are nearly Gaussian and therefore Simon's cri- 
terion is expected to be very reliable. Furthermore, at 
any t > t SI > t G the reference Gaussian state is ob- 
viously separable and thus the non-Gaussianity can be 
compared with a measure of entanglement: the relative 
entropy [24| E(g) = mm a e n[S(g\\a)] that quantifies the 
distance between g and the whole set of separable states 
f2. When t > t SI one has that 



E(g) < 6(g) « 1 



(4) 



and this confirms that in this limit the states are very 
poorly entangled (if they are) and SI allows to reliably 
estimate t s . The fact that the other criteria provide very 
close bounds on t s strengthens our conclusion. At high 
T Gaussification times are greater than all t K and we 
cannot draw the same conclusions. However, the agree- 
ment between different criteria is still an indication that 
t K may represent a good estimate of t s . 

We now focus our attention on the dependence on 
tsi on the initial non-Gaussianity 8q. Indeed, we have 
S a = 2/(eL) where f(x) = (x + 1/2) log (x + 1/2) - (x - 
1/2) log (x — 1/2) monotonically increases with x and 
cL = v/ (N + 1/2) 2 — |C| 2 is the least symplectic eigen- 
value of the covariance matrix. Upon defining g = f , 
t SI can be written as 



1, L , y/(iV + l/2) 2 -.9W2)-AM ^ 
*« = F 108^1+ Wt j (5) 

which shows that t SI is a decreasing function of Sq at any 
fixed N, and it is maximized by TWB for which 5q = 0. 
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FIG. 2: Separation times ts (straight lines) and initial non- 
Gaussianities So (dashed lines) as a function of initial energy 
N (left) and initial entanglement eo (right) for different PNES 
classes: PASV (brown, circle), PSSV (orange, square), TMC 
(red, triangle), TWB (blue) and random (black, rhombus) 
states. In both cases the bath has fixed N T = 10~ 3 . 

In the left panel of Fig. [5] the separation times and 
the initial non-Gaussianities of different PNES are plot- 
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ted against N (for t s we use t M = max K t K (g)): at any 
fixed energy N, the states with higher So have shorter 
separation times. This result holds VT and it is related 
with the fact that at any fixed N Gaussian states are 
maximally entangled Q. Let us now discuss the rela- 
tion between separation times, non-Gaussianity and the 
initial entanglement of the states. As shown in the right 
panel of Fig. [2j where *m and So are plotted as a function 
of the initial entanglement eo, the dependence is by no 
means universal. However, we notice that also at fixed 
eo states with higher Sq show shorter t s : this trend is 
not represented by an exact relation, but it represents a 
clear indication that non-Gaussianity speeds up the loss 
of entanglement, making Gaussian entanglement more 
robust than non-Gaussian one. Therefore the robustness 
of Gaussian entanglement may be conjectured to be a 
general feature of CV systems evolving in noisy Marko- 
vian channels. In fact, within the Markovian approxima- 
tion, propagation in CV noisy channels corresponds to a 
ME in Linblad form, which induces a Gaussian map and 
enforces Gaussification of any initial state. 

The results of our analysis, together with the above 
discussions, naturally lead us to formulate the following 
general conjecture: for any fixed value of the global energy 
of a PNES, and for any given noisy Markovian evolution 
with losses and thermal hopping, the Gaussian states are 
those that have maximal separation times. Besides, from 
Fig[2] we also extract another relevant feature: in the 
high-energy limit there is an approximate universality in 
separation times i.e., t s are nearly constant and simi- 
lar for all classes of states, including randomly-generated 
states, independent of the non-Gaussianity: the effect of 
the departure from Gaussianity is very small. 

Let us summarize the results of our analysis. We have 
considered a class of states (PNES) including Gaussian 
and non-Gaussian subclasses and exploited several en- 
tanglement criteria to estimate their separation times in 
a noisy channel. The analysis shows that no criterion 
is able to give better bounds than those provided by Si- 
mon's criterion. At low temperature, the estimate pro- 
vided by Simon's criterion is very reliable since PNES 
gaussify well before they lose entanglement, whereas at 
high temperature it represents a lower bound on separa- 
tion time. At any fixed energy N, separation times de- 
crease with the initial non-Gaussianity So, both at high 
and at low temperature, whereas for any fixed initial en- 
tanglement eo separation time is longer for states with 
lower So, i.e. Gaussian entanglement is the most robust 
against noise. Finally, in the high energy limit and in- 
dependently of the temperature, the differences among 
separation times of different subclasses are small, non- 
Gaussian entanglement being nearly as robust as Gaus- 
sian one. 

In conclusion, we have provided several evidences sup- 
porting the conjecture that at fixed energy Gaussian en- 



tanglement is the most robust against noise in a Marko- 
vian Gaussian channel. On the other hand, our analysis 
shows that robustness of non-Gaussian states is compa- 
rable with that of Gaussian states for sufficiently high en- 
ergy of the states. This implies that in these regimes non- 
Gaussian resources can be exploited to improve quantum 
communication protocols approximately over the same 
distances. 
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